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We measure the AjJ lifetime in the fully reconstructed decay Aj — >■ J/tphP using 10.4 fb"^ of 
pp collisions collected with the DO detector at -^i = 1.96 TeV. The lifetime of the topologically 
similar decay channel i?" — >■ J/ip Kg is also measured. We obtain t(Aj) = 1.303 ± 0.075 (stat.) ± 
0.035 (syst.) ps and t(_B") = 1.508 ±0.025 (stat.) ±0.043 (syst.) ps. Using these measurements, we 
determine the lifetime ratio of r(ASJ)/r(B'') = 0.864 ± 0.052 (stat.) ± 0.033 (syst.). 

PACS numbers: 14.20. Mr, 13.25.Hw, 13.30.Eg, 14.40. Nd 



Lifetime measurements of particles containing b quarks 
provide important tests of the significance of strong in- 
teractions between the constituent partons in the weak 
decay of b hadrons. These interactions produce mea- 
surable differences between b hadron lifetimes that the 
heavy quark expansion (HQE) [l| predicts with good ac- 
curacy through the calculation of lifetime ratios. While 
the agreement of the ratios between experimental mea- 
surements and HQE is excellent for B mesons 0, there 
are remaining discrepancies between experimental results 
and theoretical predictions for b baryons. Recently, the 
CDF Collaboration Q used the exclusive decay A° — >• 
J/%phP to report the single most precise determination 
of the h.^ lifetime which is more than 2 standard de- 
viations higher than the world average Q and slightly 
higher than the lifetime. The CDF measurement of 
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the lifetime ratio, T(A|J)/r(_B"), is higher than the HQE 
calculation including ©(l/m^) effects, 0.88 ± 0.05 [5|. 
On the other hand, theoretical predictions are in agree- 
ment with measurements by the DO Collaboration in the 
J/iIjAP and semileptonic 0] channels, by the CDF Col- 
laboration in the A+tt" final state Q, by the DELPHI, 
OPAL and ALEPH Collaborations in semileptonic de- 
cays 0-111, and previous measurements also in semilep- 
tonic channels by the CDF Collaboration • More mea- 
surements of the A^ lifetime and of the ratio r(A°)/T(i?°) 
are required to resolve this discrepancy. 

In this article we report a measurement of the A|j life- 
time using the exclusive decay A^ — J /ip A". The i?° life- 
time is also measured in the topologically similar chan- 
nel B^ — J/ipKg. This provides a cross-check of the 
measurement procedure, and allows the lifetime ratio to 
be determined directly. The data used in this analysis 
were collected with the DO detector during the complete 
Run II of the Tevatron Collider, from 2002 to 2011, and 
correspond to an integrated luminosity of 10.4 fb"^ of 
pp collisions at a center of mass energy y/s = 1.96 TeV. 

A detailed description of the DO detector can be found 
in Refs. [13-16j. Here, we describe briefly the most rel- 
evant detector components used in this analysis. The 
DO central tracking system is composed of a silicon mi- 
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crostrip tracker (SMT) and a central scintillating fiber 
tracker (CFT) immersed in a 2 T solenoidal field. The 
SMT and the CFT are optimized for tracking and vertex- 
ing for the pseudorapidity region \r]\ < 3.0 and \ri\ < 2.0, 
respectively, where rj = — ln[tan(6'/2)] and 9 is the po- 
lar angle with respect to the proton beam direction. 
Preshower detectors and electromagnetic and hadronic 
calorimeters surround the tracker. A muon spectrometer 
is located beyond the calorimeter, and consists of three 
layers of drift tubes and scintillation trigger counters cov- 
ering |?7| < 2.0. A 1.8 T toroidal iron magnet is located 
outside the innermost layer of the muon detector. 

For all Monte Carlo (MC) simulations in this arti- 
cle, we use PYTHIA [l3| to simulate the pp collisions, 
EVTGEN [3] for modeling the decay of particles contain- 
ing b and c quarks, and GEANT [I^ to model the detector 
response. Multiple pp interactions are modeled by over- 
laying hits from random bunch crossings onto the MC. 

In order to reconstruct the A° and candidates, we 
start by searching for J/ip — )■ /x^/x" candidates, which 
are collected by single muon and dimuon triggers. The 
triggers used do not rely on the displacement of tracks 
from the interaction point. At least one pp interaction 
vertex (PV) must be identified in each event. The in- 
teraction vertices are found by minimizing a function 
that depends on all reconstructed tracks in the event and 
uses the transverse beam position averaged over multiple 
beam crossings. The resolution of the PV is ~ 20 fiui in 
the plane perpendicular to the beam (transverse plane). 
Muon candidates are reconstructed from tracks formed 
by hits in the central tracking system and with trans- 
verse momentum {px) greater than 1 GeV/c. At least 
one muon candidate in the event must have hits in the 
inner layer, and in at least one outer layer of the muon de- 
tector. A second muon candidate, with opposite charge, 
must either be detected in the innermost layer of the 
muon system or have a calorimeter energy deposit con- 
sistent with that of a minimum-ionizing particle along 
the direction of hits extrapolated from the central track- 
ing system. Each muon track is required to have at least 
2 hits in the SMT and 2 hits in the CFT to ensure a high 
quality common vertex. The probability associated with 
the vertex fit must exceed 1%. The dimuon invariant 
mass is required to be in the range 2.80 — 3.35 GeV/c^, 
consistent with the J/ip mass. 

Events with J/ip candidates are reprocessed with a ver- 
sion of the track reconstruction algorithm that increases 
the efficiency for tracks with low pt and high impact 
parameter [20|. We then search for A° — >■ pTr~ can- 
didates reconstructed from pairs of oppositely charged 
tracks. The tracks must form a vertex with a probabil- 
ity associated with the vertex fit greater than 1%. The 
transverse impact parameter significance (the transverse 
impact parameter with respect to the PV divided by its 
uncertainty) for the two tracks forming A° candidates 
must exceed 2, and 4 for at least one of them. Each 



A*^ candidate is required to have a mass in the range 
1.105 - 1.127 GeV/c^ The track with the higher pr is 
assigned the proton mass. MC simulations indicate that 
this is always the correct assumption, given the track 
Pt detection threshold of 120 MeV/c. To suppress con- 
tamination from decays of more massive baryons such 
as S° — > A''7 and S° — >■ A°7r°, the A° momentum vec- 
tor must point within 1 degree back to the J/ip vertex. 
The same selection criteria are applied in the selection 



of K' 



candidates, except that the mass win- 



dow is chosen in the range 0.470 — 0.525 GeV/c^ and 
pion mass assignments are used. Track pairs simultane- 
ously reconstructed as both A° and Kg, due to different 
mass assignments to the same tracks, are discarded from 
both samples. This requirement rejects 23% (6%) of the 
A^ -> J/tP A° — >• J/ip K^) signal, as estimated from 
MC, without introducing biases in the lifetime measure- 
ment. The fraction of background rejected by this re- 
quirement is 58% (48%) as estimated from data. It is 
important to remove these backgrounds from the sam- 
ples to avoid the introduction of biases in the lifetime 
measurements. 

The A° candidates are reconstructed by performing a 
kinematic fit that constrains the dimuon invariant mass 
to the world-average J/ip mass [3], and the A° and two 
muon tracks to a common vertex, where the A*' has been 
extrapolated from its decay vertex according to the re- 
constructed A'^ momentum vector. The invariant mass 
of the A° candidate is required to be within the range 
5.15 — 6.05 GeV/c^. The PV is recalculated excluding the 
A^ final decay products. The final selection requirements 
arc obtained by maximizing S = S/ ^ S + B, where S 
{B) is the number of signal (background) candidates in 
the data sample: the decay length of the A" (measured 
from the A|J vertex) and its significance are required to be 
greater than 0.3 cm and 3.5, respectively; the pt of the 
J A° and A*^ daughter tracks are required to be greater 
4.5, 1.8 and 0.3 GeV/c, respectively; and the isolation of 
the A° [2lj is required to be greater than 0.35. After this 
optimization, if more than one candidate is found in the 
event, which happens in less than 0.3% of the selected 
events, the candidate with the best A° decay vertex fit 
probability is chosen. We have verified that this selection 
is unbiased by varying the selection values chosen by the 
optimization as described in more detail later. The same 
selection criteria are applied to i?" — ?> J/iJ^Kg decays, 
except that the B'^ mass window is chosen in the range 
4.9-5.7 GeV/c2. 

The samples of A[j and B^ candidates have two pri- 
mary background contributions: combinatorial back- 
ground and partially reconstructed b hadron decays. The 
combinatorial background can be divided in two cate- 
gories: prompt background, which accounts for « 70% 
of the total background, primarily due to direct produc- 
tion of J/ip mesons; and non-prompt background, mainly 
produced by random combinations of a Jf-ip meson from 
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FIG. 1: (color online) Invariant mass distributions for (a) AjJ — >■ J/iplSP and (b) £?" — >■ J/ipRg candidates, with fit results 
superimposed. Events in mass regions contaminated with partially reconstructed b hadrons (hatched region) are excluded from 
the maximum likelihood function used to determine the and B° lifetimes. 



a b hadron and a A° (Kg) candidate in the event. Con- 
tamination from partially reconstructed b hadrons come 
from b baryons (B mesons) decaying to a J/ip meson, a 
A° baryon {Kg meson), and additional decay products 
that are not reconstructed. 

We define the transverse proper decay length as A = 
cMLxy/pT, where M is the mass of the b hadron taken 
from the PDG [4], and L^y is the vector pointing from 
the PV to the b hadron decay vertex projected on the 
b hadron transverse momentum (px) direction. Due to 
the fact that signal and partially reconstructed b hadron 
decays have similar A distributions that are particularly 
hard to disentangle in the lifetime fit, we remove par- 
tially reconstructed b hadrons by rejecting events with 
Ag (5°) invariant mass below 5.42 (5.20) GeV/c^ from 
the A° (B^) sample, as shown in Fig. [TJ This figure 
shows the and B^ invariant mass distributions with 
results of unbinned maximum likelihood fits superim- 
posed, excluding events in zones contaminated by par- 
tially reconstructed b hadrons. The signal peak is mod- 
eled by a Gaussian function. The combinatorial back- 
ground is parametrized by an exponentially decaying 
function, while partially reconstructed b hadrons are de- 
rived from MC. It can be seen from Fig. [1] that partially 
reconstructed b hadrons contribute minimally to the sig- 
nal mass region. 

In order to extract the lifetimes, we perform separate 
unbinned maximum likelihood fits for AJJ and candi- 
dates. The likelihood function (£) depends on the prob- 
ability of reconstructing each candidate event j in the 
sample with the mass rrij, the proper decay length Xj 
and proper decay length uncertainty aj*: 

>C =n [.fsJ's{m,,X„<j^) + {l~fs)Mm„ X^'jf)] , (1) 



where fs is the fraction of signal events, and J-g {Fb) is 
the product of the probability distribution functions that 
model each of the three observables being considered for 
signal (background) events. The background is further 
divided into prompt and non-prompt components. For 
the signal, the mass distribution is modeled by a Gaus- 
sian function; the A distribution is parametrized by an 
exponential decay, e^'^^/^'^/cr, convoluted with a Gaus- 
sian function TZ = er '^^'"^'i / \/2t: sa^ that models the 
detector resolution; the distribution is obtained from 
MC simulation and parametrized by a superposition of 
Gaussian functions. Here r is the lifetime of the b hadron, 
and the event- by- event uncertainty is scaled by a 
global factor s to take into account a possible underesti- 
mation of the uncertainty. For the background, the mass 
distribution of the prompt (non-prompt) component is 
modeled by a constant (exponential) function as observed 
in data when the requirement A > 100 /im is imposed; the 
prompt component of the A distribution is parametrized 
by the resolution function, and the non-prompt compo- 
nent by the superposition of two exponential decays for 
A < and two exponential decays for A > 0, as observed 
from events in the high-mass sideband of the b hadron 
peak (above 5.80 and 5.45 GeV/c^ for A|] and B° , re- 
spectively). Finally, the background distribution is 
modeled by two exponential functions convoluted with 
a Gaussian function as determined empirically from the 
high-mass sideband region. In total, there are 19 pa- 
rameters in each likelihood fit: lifetime, mean and width 
of the signal mass, signal fraction, prompt background 
fraction, one non-prompt background mass parameter, 7 
non-prompt background A parameters, 5 background 
parameters, and one resolution scale factor. 

The maximum likelihood fits to the data yield 
cr(Ag) = 390.7 ± 22.4 /xm and cr(B°) = 452.2 ± 7.6 ^jluy. 
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The numbers of signal events, derived from are 
755 ± 49 (Ag) and 5671 ± 126 Figure [2] shows the 

A distributions for the A° and the _B" candidates. Fit 
results are superimposed. 

We investigate possible sources of systematic uncer- 
tainties on the measured lifetimes related to the models 
used to describe the mass, A, and cr^ distributions. For 
the mass we consider a double Gaussian to model the 
signal peak instead of the nominal single Gaussian, an 
exponential function for the prompt background in place 
of a constant function, and a second order polynomial 
for the non-prompt background. The alternative mass 
models are combined in a single maximum likelihood fit 
to take into account correlations between the effects of 
the different models, and the difference with respect to 
the result of the nominal fit is quoted as the systematic 
uncertainty on the mass model. For A we study the fol- 
lowing variations: the introduction of a second Gaussian 
function along with a second scale factor to model the 
resolution, the exponential functions in the non-prompt 
background replaced by exponentials convoluted with the 
resolution function, one non-prompt negative exponen- 
tial instead of two, and one long positive exponential 
together with a double-Gaussian resolution as a substi- 
tute for two non-prompt exponentials and one Gaussian 
resolution. All A model changes are combined in a fit, 
and the difference between the results of this fit and the 
nominal fit is quoted as the systematic uncertainty due 
to A parametrization. For we use two different ap- 
proaches: we use the distribution extracted from data by 
background subtraction, parameterized similarly to the 
nominal background a"** model, instead of the MC model, 
and we use cr"^ distributions from MC samples generated 
with different {B^) lifetimes. The largest variation in 
the lifetime (with respect to the nominal measurement) 
between these two alternative approaches is quoted as 
the systematic uncertainty due to cr^ parametrization. 
Residual effects due to contamination from partially re- 
constructed b hadrons in the samples are investigated by 
changing the requirement on the invariant mass of the 
AJJ and i?" candidates which are included in the likeli- 
hood fits: the threshold is moved to lower (higher) in- 
variant masses by 40 (20) MeV/c^, where 40 MeV/c^ is 
the resolution on the invariant mass of the reconstructed 
signal. The largest variation in the lifetime is quoted as 
the systematic uncertainty due to possible contamination 
from partially reconstructed b hadrons. In the lifetime 
fit the contamination from the fully reconstructed decay 
— )• J/ipKg is assumed to have little impact on the 
final result. To test this assumption the B'^ J/ipKg 
contribution is included in the non-prompt component. 
The lifetime shift is found to be negligible. The sys- 
tematic uncertainty due to the alignment of the SMT 
detector was estimated in a previous study 0] by recon- 
structing the i?° sample with the positions of the SMT 
sensors shifted outwards radially by the alignment uncer- 



TABLE I: Summary of systematic uncertainties on the mea- 
surements of cr(Aj) and cr{B°), and on their ratio. Individ- 
ual uncertainties are combined in quadrature to obtain the 
total uncertainties. 



Source 


AO (/im) 


B° (^m) 


Ratio 


Mass model 


2.2 


6.4 


0.008 


Proper decay length model 


7.8 


3.7 


0.024 


Proper decay length uncertainty 


2.5 


8.9 


0.020 


Partially reconstructed b hadrons 


2.7 


1.3 


0.008 






0.4 


0.001 


Alignment 


5.4 


5.4 


0.002 


Total 


10.4 


12.9 


0.033 



tainty and then fitting for the lifetime. The systematic 
uncertainties are summarized in Table HI 

We perform several cross-checks of the lifetime mea- 
surements. We extract the signal yield in bins of A by 
fitting the mass distribution in each of these regions. 
From these measurements, lifetimes are obtained by the 
minimization of the signal yield expected in each A 
bin according to the first term in Eq. [TJ While this 
method is statistically inferior with respect to the maxi- 
mum likelihood fit, it is also less dependent on the mod- 
eling of the different background components. The re- 
sults of this study are cr^o = 391.4 ±35.8 (stat.) /im and 
ctbo = 458.3 ± 8.9 (stat.) /im. The sample is also split 
into different data taking periods, rj regions, and num- 
bers of hits in the SMT detector. All results obtained 
with these variations are consistent with our measure- 
ment. In order to check that the optimization procedure 
does not give a potential bias to the selection, we verify 
that our results remain stable when all requirements in 
variables used in the optimization process are removed 
one at a time, when looser and tighter requirements are 
applied to kinematic variables, and when multiple candi- 
dates that pass all selection requirements per event are 
allowed. The results also remain stable after removing 
the high-end tail (above 100 /im) of the distribution, 
mainly populated by background events. We also cross 
check the fitting procedure and selection criteria by mea- 
suring the A" and 5° lifetimes in MC events. The life- 
times obtained are consistent with the input values. 

In summary, using the full data sample collected by the 
DO experiment, we measure the lifetime of the Aj] baryon 
in the J /ipA^ final state to be 

T(Ag) = 1.303 ± 0.075 (stat.) ± 0.035 (syst.) ps, (2) 

consistent with the world-average, 1.425 ± 0.032 ps 0]. 
The method to measure the A° lifetime is also used for 
B^ — > J/tpKg decays, for which we obtain 

t{B°) = 1.508 ± 0.025 (stat.) ± 0.043 (syst.) ps, (3) 

in good agreement with the world average, 1.519 ± 
0.007 ps (1. 
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(a) D0, 10.4fb"' (b) D0, 10.4fb"' 




FIG. 2: (color online) Proper decay length distributions for (a) — >■ J/ip hP and (b) B'^ — > J/tpKg candidates, with fit results 
superimposed. 



Using these measurements we calculate the ratio of life- 
times, 

t(A°) 

-^4t = 0-864 ± 0.052 (stat.) ± 0.033 (syst.), (4) 

where the systematic uncertainty is determined from the 
differences between the lifetime ratio obtained for each 
systematic variation and the ratio of the nominal mea- 
surements, and combining theses differences in quadra- 
ture, as shown in Table ID Our result, 0.86 ± 0.06, is in 
good agreement with the HQE prediction of 0.88±0.05 l5\ 
and compatible with the current world-average, 1.00 ± 
0.06 1^], but differs with the latest measurement of the 
CDF Collaboration, 1.02 ±0.03 >], at the 2.2 standard 
deviations level. Our measurements supersede the previ- 
ous DO results of r(Ag), t(B") and r(Ag)/T(B°) ;6]. 
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